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Abstract — Second order statistics provides a dynamic repre- 
sentation of a fading channel and plays an important role in the 
evaluation and design of the wireless communication systems. 
In this paper, we present a novel analytical framework for the 
evaluation of important second order statistical parameters, as 
the level crossing rate (LCR) and the average fade duration 
(AFD) of the amplify-and-forward multihop Rayleigh fading 
channel. More specifically, motivated by the fact that this channel 
is a cascaded one and can be modeled as the product of N fading 
amplitudes, we derive novel analytical expressions for the average 
LCR and the AFD of the product of N Rayleigh fading envelopes 
(or of the recently so-called iV*Rayleigh channel). Furthermore, 
we derive simple and efficient closed-form approximations to the 
aforementioned parameters, using the multivariate Laplace ap- 
proximation theorem. It is shown that our general results reduce 
to the corresponding ones of the specific dual-hop case, previously 
published. Numerical and computer simulation examples verify 
the accuracy of the presented mathematical analysis and show 
the tightness of the proposed approximations. 

Index Terms — Multihop relay communications, level crossing 
rate, average fade duration, Laplace approximation, Rayleigh 
fading. 

I. Introduction 

MULTIHOP communications have recently emerged as 
a viable option for providing broader and more eff- 
icient coverage both in traditional (e.g. bent pipe satellites) 
and modern (e.g. ad-hoc, WLAN) communications networks. 
In such systems, contrary to conventional wireless networks, 
several intermediate terminals operate as relays between the 
source and the destination [ 1 ]-[ 1 3] . 

Multihop transmissions can be categorized as either 
non-regenerative (amplify-and-forward, AF) or regenerative 
(decode-and-forward, DF), depending on the relay function- 
ality. In DF systems, each relay decodes its received signal 
and then re-transmits this decoded version. On the other hand, 
in AF systems, the relays just amplify and re-transmit their 
received signal. Furthermore, the AF systems can use either 
channel state information (CSI)-assisted relays [1] or fixed- 
gain relays [2] (also known as blind or semi-blind relays 
[10]). A (CSI)-assisted relay uses instantaneous CSI of the 
channel between the transmitting terminal and the receiving 
relay terminal to adjust its gain, whereas a fixed-gain relay just 
amplifies its received signal by a fixed gain [2] [10]. Note, that 
systems with fixed-gain relays perform close to the systems 
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with (CSI)-assisted relays [2], while their easy deployment and 
low complexity make them attractive from a practical point of 
view. 

Several works in the open literature have provided perfor- 
mance analysis of AF and DF systems in terms of bit error rate 
(BER) and outage probability under different assumptions for 
the amplifier gain [1]-[13]. Among them, only two works deal 
with the dynamic, time-varying nature of the underlying fading 
channel [12], [13], despite the fact that it is necessary for the 
system's design or rigorous testing [ 1 3]-[ 15] . In [12], the level 
crossing rate (LCR) and the average fade duration (AFD) of 
multihop DF communication systems over generalized fading 
channels is studied, both for noise-limited and interference- 
limited systems, while Patel et. al in [13] provide useful 
exact analytical expressions for the AF channel's temporal 
statistical parameters, such as the auto-correlation and the 
LCR. However, the approach presented in [13] is limited only 
to the dual-hop fixed-gain AF Rayleigh fading channel. 

In this paper, we study the second order statistics of the 
multihop fixed-gain AF Rayleigh fading channel. More specif- 
ically, motivated by the fact that this channel is a cascaded one 
and can be modelled as the product of N fading amplitudes, 
we derive a novel analytical framework for evaluation of the 
average LCR and the AFD of the product of N Rayleigh 
fading envelopes. Since the presented exact expressions are 
computationally attractive only for small values of N, we 
derive simple and yet efficient closed-form approximations 
using the multivariate Laplace approximation theorem [22, 
Chapter IX. 5], [23], which can be efficiently used to evaluate 
the aforementioned second order statistical parameters. These 
important theoretical results are then applied to investigate 
the second order statistics of the multihop Rayleigh fading 
channel. Numerical and computer simulations validate the 
accuracy of the presented mathematical analysis and show the 
tightness of the proposed approximations. 

The remainder of the paper is organized as follows: In the 
next section, the second order statistics analysis of the product 
of N Rayleigh fading amplitudes is presented, providing exact 
and tight approximated expressions for the LCR and the AFD. 
In section IV, these theoretical results are applied to study the 
LCR and the AFD of the fixed-gain multihop relay fading 
channel. Section V includes numerical and simulations results 
which validate the proposed mathematical analysis while some 
concluding remarks are given in section VI. 

II. Level Crossing Rate and Average Fade 
Duration of The Product of N Rayleigh Envelopes 

Let {Xi(t)}fL 1 be N independent and not necessarily 
identically distributed (i.n.i.d.) Rayleigh random processes, 
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each distributed according to [16], [17], 

,.2 ' 



JXi [x) = — exp 



x 



x > 0, 



(1) 



in an arbitrary moment t, where Q, — E{Xf(t)} is the mean 
power of the i-th random process (1 < i < N), with E{-} 
denoting expectation. 

If {Xi(t)}fL 1 represent received signal envelopes in an 
isotropic scattering radio channel exposed to the Doppler 
effect, they must be considered as time-correlated random 
processes with some resulting Doppler spectrum. This Doppler 
spectrum differs depending on whether a fixed-to-mobile chan- 
nel [16]-[17] or mobile-to-mobile channel [ 1 8]-[ 19] appears 
in the particular wireless communication system. However, 
in both cases, it was found that the time derivative of i-th 
envelope is independent from the envelope itself, and follows 
the Gaussian probability distribution function (PDF) [16]-[19] 
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with variance 
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If envelope Xi is formed on a fixed-to-mobile channel, then 
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where f m i is the maximum Doppler frequency shift induced 
by the motion of the mobile station [16]-[17]. If envelope Xi 
is formed on a mobile-to-mobile channel, then 
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where f mi and f mi are the maximum Doppler frequency 
shifts induced by the motion of both mobile stations (i.e., the 
transmitting and the receiving stations, respectively) [19]. It is 
important to underline that the maximum Doppler frequency 
in a fixed-to-mobile channel is 



fd max f n 



(6) 



whereas the maximum Doppler frequency in a mobile-to- 
mobile channel is 



fdi 



f mi f rr 
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The above results are essential in deriving the second-order 
statistical parameters of individual envelopes, as the average 
LCR and the AFD [16], [17], [19]. 

Below, we derive exact analytical and approximate solutions 
for both of the above parameters for the product of the N 
Rayleigh envelopes, 



Y(t)=l[X t (t), 



(8) 



where Y(t) in $Q is iV*Rayleigh random process or, at any 
given moment t, iV*Rayleigh random variable, following the 
definition given in [20]. 



For specified values {Xi}fL 1 — {a^}^, the product Y is 
fixed to the specific value 
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The LCR of Y at threshold y is defined as the rate at which 
the random process crosses level y in the negative direction 
[16]. To extract LCR, we need to determine the joint PDF 
between Y and Y, fyyiVi y)> an d to apply the Rice's formula 
[17, Eq. (2.106)], 
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Our method does not require explicit determination of 
/yy(j/,y) in order to obtain analytically the LCR of the 
./V*Rayleigh random process, as presented below. 

First, we need to find the time derivative of dH), which is 
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Conditioning on the first N — 1 envelopes {Xi\^ =l l = 
{xi}^ 1 , we have the conditional joint PDF Y and Y written 

as fYY\x 1 ---x N - 1 (y'y\ Xl '---> XN - i y This conditional joint 
PDF can be averaged with respect to the joint PDF of the 
N — 1 envelopes {^i}^ 1 t0 produce the required joint PDF 
of Y and Y as 

fyviy^y) 
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where to derive ( TTZb the mutual independence of the N — 1 
envelopes is used. 

The conditional joint PDF / y y| Xl ... Xjv _ 1 (y, y\xi, 
xjv-i) can be further simplified by fixing Y = y and using 
the total probability theorem, 



/YYIAV-.A-jv-^y^Nl' -,XN-l) 

= fY\YX 1 -X N - 1 (v\V,H!l,...,X N - 1 ) 

x /y|Xi-JfAr_i(y|a!i,- 



,XiV-i), (13) 



where each of the two multipliers in dT3l can be determined 
from the above defined individual PDFs and their parameters. 

Based on ( fTTT i. the conditional PDF fyiYXi—XN livlVi Xl > 
...,xn-i) is easily established to follow the Gaussian PDF 
with zero mean and variance 



Y\YX 




(14) 



3 



The conditional PDF of Y, given {Xi}^ 1 = {xi}^ 1 , 
that appears in ( fT~3b is easily determined in terms of the PDF 
of the remaining A^-th envelope, 
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Introducing (\3$ and (|T5j into (|T2j, then (O into dTOJ, and 
changing the orders of the integration, we obtain 
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The bracketed integral in ( TTBT l is found using (fl4l > as 
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By substituting (HJ and ( TTTb into ( TToT ). we obtain the final 
exact formula for the LCR as 
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In principle, JT8l together with JT9] l provide an exact 
analytical expression for the LCR of the product of the product 
of N Rayleigh envelopes (i.e., 7V*Rayleigh random process 
[20]). However, ( TT~8T > becomes computationally attractive only 
for small values of N, such as N — 2 and N = 3, where it is 
possible to apply multidimensional numerical integration (as 
Gaussian-Hermite quadrature [27]), included in most of the 
well-known mathematical software packages. 

The AFD of Y at threshold y is defined as the average time 
that the 7V*Rayleigh random process remains below level y 
after crossing that level in the downward direction, 
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where Fy(-) denotes the cumulative distribution function 
(CDF) of Y. Fortunately, Fy (•) was derived recently in closed- 
form [20, Eq. (7)], as 



F Y (y) = G$ +1 
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where G[] is the Meijer's G-function [21, Eq. (9.301)]. Note 
that Meijers G-function is a standard built-in function in well- 
known mathematical software packages, such as MAPLE and 
MATHEMATICA. 

A. An Approximate Solution for the LCR 

Next, we present a tight closed-form approximation of ( TT~8b 
using the multivariate Laplace approximation theorem [22, 
Chapter IX. 5], [23] for the Laplace-type integral 



J(A) = 



u(x) exp(—Xh(x))dx, 



(22) 
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where u and h are real-valued multivariate functions of x = 
[xi, ■ ■ ■ , xn-i], A is a real parameter and D is unbounded 
domain in the multidimensional space 
A comparison of (fT8l and d22l yields 



u(x) = 



h(x) 




(23) 



(24) 



and A = 1. A brief description of the multivariate Laplace 
approximation theorem and its applicability conditions are 
provided in the Appendix A. 

Note, that in the case of ( fl~8l ), all the applicability conditions 
of the theorem are fulfilled. Namely, within the domain of 
interest D, the function h(x) has a single interior critical point 
x = [x\, ■ ■ ■ jiGjv-i], where 
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which is obtained from solving the set of equations dh/dx,i = 
0, where 1 < i < N - 1. The Hessian (N - 1) X (N - 1) 
square matrix A, defined by (A.2), is written as 



8/fti 



fiiQ 
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f^^TV-l 

8/Qn-i 

(26) 

By using induction, it is easy to determine that the N — 1 
eigenvalues of A are calculated as fa = 4/0^ for 1 < i < 
N — 2, and ^jv_i = 4AT/f2jv-i- Thus, all eigenvalues of A 
are positive, which, by definition, means that the matrix A 
is positive definite. By means of the second derivative test, 
since the Hessian matrix A is positive definite at point x, 
h{x) attains a local minimum at this point (which in this case 
is the absolute minimum in the entire domain D). 
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At this interior critical point x, it holds 
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where (l27l i is obtained using (0). Thus, by using (A. 3), it is 
possible to approximate d22b for large A as 



(iV-l)/2 




jB. Special Cases 

If the mean powers of all envelopes are assumed mutually 
equal, {f^}^ = f2, (l3TT l reduces to 

1/2 
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which is the approximate closed-form solution for the LCR of 
the product of N identically distributed Rayleigh envelopes. 

Interestingly, when N = 1, rfiil ) further reduces to the 
classic expression for the LCR of a Rayleigh-faded signal, 
regardless of its mean power f2, i.e. 



Ny{y) = h ^j= ^[- V - 



(34) 



where f\ = f m for the fixed-to-mobile channel [16]- [17], and 
( - 29 - ) fx = V 1m + fm for the mobile-to-mobile channel [19]. 



It is well-know that the determinant of the square matrix is 
equal to the product of its eigenvalues, so det(A) is calculated 

as 

7Y2 2 ( JV - 1 ) fi Jv iV2 2 ( iV - 1 ) 



det(A) 
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Although approximation ( 1291 is proven for large A [22] -[23], 
it is often applied when A is small and is observed to be very 
accurate as well. Similarly to [24], we apply the theorem for 
A = 1. Therefore, the approximate closed-form solution for 
the LCR of jV*Rayleigh random process Y at threshold y is 
determined by 
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The numerical results presented in Section IV validate the 
high accuracy of the Laplace approximation applied for our 
particular case. 

Combining (ED and CD) into the AFD of N Rayleigh 
random process Y at threshold y is approximated as 
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III. Second Order Statiscs Of Multihop 
Transmissions 

Next, we apply the theoretical results of the previous Section 
to analyze the second order statistics of the multihop Rayleigh 
fading channel. 

A. System Model 

Let's consider a multihop wireless communications system, 
operating over i.n.i.d. flat fading channels (Fig. 1). The source 
station S communicates with the destination station D through 
jV — 1 nodes Ti, T%, T^-i, which act as intermediate 
relays from one hop to the next. These intermediate nodes 
are employed with non-regenerative relays with fixed gain Gi 
given by 



G = 



1 



(35) 



with Go = 1 and Cq = 1 for the source S. In d35l) . Wo,i is 
the variance of the Additive White Gaussian Noise (AWGN) 
at the output of the i-th relay, and d is a constant for the 
fixed-gain Gi. 

Assume that terminal S is transmitting a signal s(t) with an 
average power normalized to unity. Then, the received signal 
at the first intermediate node, T%, at moment t, can be written 

as 

rx{t) = a 1 (t)s{t)+wi(t), (36) 

where a\{t) is the fading gain between S and T\, and w\{t) 
is the AWGN at the input of T\ with variance Wq,x- The 
signal r\ is then multiplied by the gain G\ of the node T\ 
and re-transmitted to node T2, where its received signal can 
be written as 



r a (f) - Gxa 2 (t) (ai (t)s(t) + wx{t)) + w 2 {t) 



(37) 



with 0.2(f) being the fading gain of the channel between T\ 
and T2. Generally, the received signal at the fc-th relay Tk 
(k = 1, 2, jV - 1) is given by 



r k (t) = G k ~xctk(t)r k -i(t) + w k (t) 



(38) 
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Fig. 1. Multihop wireless communication system 



finally resulting in a total fading gain at the destination station 
D, given by 



N 



a(t) = Y[ ai (t)G 



i-l 



(39) 



B. LCR and AFD of Multihop Transmissions 

If the fading amplitude received at node Ti, ai(t), is a 
time-correlated (due to mobility of Ti_i and/or Ti) Rayleigh 
random process, distributed according to (HJ at any given 
moment t, with mean power 



n^E{a 2 (t)}, 
then the i-th element of the product in 

Xi(t) = ai (t)Gi-i 



(40) 



(41) 



is again a time-correlated Rayleigh random process, distributed 
according to ((TJ, with mean power 



(42) 



Comparing (JHJ and (1391 , we realize that the total fading 
amplitude at the destination station D (i.e., the received desired 
signal without the AWGN) is described as the 7V*Rayleigh 
random process Y(t) = a(t), whose average LCR and AFD 
are determined in the previous Section. 

Based on the system model from Fig. 1, if all stations 
are assumed mobile with maximum Doppler frequency shifts 
fmS, fmD , /mi(l < i < N — 1) for source S, destination D 
and relays, respectively, then for the i-th hop 



fi f m(i— 1) fmi 

with f m0 = f m s and f mN = f mD , and 

N N-l 
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Using (1181 . (l20b and (|2TT i. it is now straightforward to obtain 
the exact expressions for the average LCR and the AFD of the 



multihop Rayleigh fading channel, as follows 
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respectively, where ff is given by (T43b and $ is given by ( fT9l . 

It must be noted here that, for N = 2, d45l > is transformed 
into a single integral of the form 



N a (a) 



^2/2 +^1/1^4 exp 



a 2 x 2 
+ fh 



dx, (47) 



which, after changing integration variable x with new variable 
t according x = a ft, reduces to the known result [13, Eq. 
(17)]. 

By combining (l44l with d3~TT ) and ( |44T > with d32l . we also 
obtain approximate solutions for the average LCR and AFD 
of the multihop Rayleigh fading channel as 



N a (a) 



and 
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fmS + ^ £ fmi + /m_D 
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(271)^/2 
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a / a 2 l N 
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2 J]] /mi + /mD J 
i=l J . 




a VN 

exp I N- 



$l/N ^ 

(49) 



respectively, where $ is given by (fT9l ). 

We see that d48l and j49l approximate the average LCR and 
AFD of the total fading amplitude for arbitrary power of the 
fading amplitudes Cli, arbitrary relay gains G; and arbitrary 
maximal Doppler shifts for the nodes / m j. 

It must be noted here that, for N = 2, d48l is an efficient 
closed-form alternative to the corresponding one given by 
[13, Eq. (17)] for the dual-hop case. Furthermore, as it will 
be shown in the next section, the proposed approximation is 
highly accurate. 
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C. Special Cases 

If we assume that i) all stations are mobile and induce same 
maximal Doppler shifts (i.e., f mS = f mi = f mD = f m ), 
ii) the fading amplitudes in all hops have equal powers (i.e., 
f2j = fl), then (l48l i reduces to 



N a (a) » 72/V 



(27T) jV / 2 
' $1/2 



exp 



-N 



a 2 / N 



(50) 



where, according to CLUl, $ = 0^ I]i=i G ?- 

If we assume that i) the destination station is fixed (i.e., 
fmD = 0) but all other stations are mobile inducing same 
maximal Doppler shifts (i.e., f m s — fmi — fm), and ii) the 
fading amplitudes in all hops have equal powers (i.e., fj, = SI), 
then (l48l l reduces to 



N a (a) « /„ 



2A-1\ 1/2 (271)^/2 



A 



$1/2 



■ exp 



-iV 



(51) 



where again $ = Sl w ]X=7 G 



IV. Numerical Results and Discussion 

In this Section, we provide some illustrative examples for 
the average LCR and AFD of the fading gain process of 
the received desired signal at the destination of a multihop 
non-regenerative relay transmission system from Fig. 1. The 
numeric examples obtained from the derived approximate 
solutions are validated by extensive Monte-Carlo simulations 
over the system model described in Section III . 

Based on the system model from Fig. 1, we considered a 
multihop transmission system consisted of a source terminal 
S, 4 fixed-gain relays, and a destination terminal D. The 
destination D has fixed position, whereas the source and all 
relays are mobile and induce same maximal Doppler shift f m . 

The fixed-gain relays are assumed semi-blind with gains in 
Rayleigh fading channel calculated according to [2, Eq. (15)] 
and [11, Eq. (19)] 



G 



i.sb 



l 



exp 



r o 



(52) 



where ji — Cli/Woj is the mean SNR on the i-th hop, 
and r(-,-) is the incomplete Gamma function. Relay gain 
calculated according to ( |52l assures mean power consumption 
equal to that of a CSI-assisted relay, whose gain inverts the 
fading effect of the previous hop while limiting the output 
power at moments with deep fading. 

Depending on the stations' mobility, we used two different 
2D isotropic scattering models for the Rayleigh radio channel 
on each hop of the multihop transmission system. For the 
fixed-to-mobile channel (hop), we used the classic Jakes 
channel model [16]-[17]. For the mobile-to-mobile channel 
(hop), we used the Akki and Habber's channel model [18]- 
[19]. The Monte-Carlo simulations of the latter were realized 
by using the sum-of-sinusoids method proposed in [25] -[26]. 

Each presented figure depicts the received signal's normal- 
ized LCR (N a / fm) or normalized AFD (T a f m ) versus the 
normalized threshold (a/y/fl) at 3 different nodes along the 
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Fig. 2. Average level crossing rate, when 71 = 72 = 73 = 74 = 75 
7 = 5 dB 
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Fig. 3. Average fade duration, when 71 = 72 = 73 = 74 = 75 = 7 = 5 
dB 



multihop transmission system: at relay T2 (curve denoted by 
N = 2), at relay T3 (curve denoted by = 3) and at the 
destination D (curve denoted by N = 5). Note that, when 
applying the considered scenarios in d48l > or dBTT l. a and fl 
appear together as a/ y/Cl. 

Figs. 2-5 assume equal power of the fading amplitudes in all 
hops fti = SI, and equal variance of the AWGN Wo t i = Wq. 
Thus, 7.; = 7, Gi.sb = G s b for 1 < i < 5, so the mean of 
Rayleigh random process Xi(t) — ai(t)G,_i !S 6 is calculated 
as 

Cli = exp (J-^J T ^0, = SI, 2 < i < 5, (53) 

whereas fii = SI is selected independently from the AWGN, 
since Gq — 1. In this case, 



$ = SI exp 



N-l 

7 



r j 0, i 
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-1 N-l 



(54) 



Figs. 6-7 assume equal powers of the fading amplitudes 
in all hops S% = SI, and unequal variances of the AWGN 
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Fig. 4. Average level crossing rate, when 71 
7 = 20 dB 



72 = 73 = 74 = 75 = Fig. 6. Average level crossing rate, when 71 = dB, 72 = 10 dB, 73 = 15 
dB, 74 = 15 dB, 75 = 20 dB 




simulation 
approximation! 




normalized threshold 



normalized threshold 



Fig. 5. 
dB 



Average fade duration, when 71 = 72 = 73 = 74 = 75 = 7 = 20 



Fig. 7. Average fade duration, when 71 
74 = 15 dB, 75 = 20 dB 



dB, 72 = 10 dB, 73 = 15 dB, 



Wo t i. Thus, the mean of Rayleigh random process Xi(t) = 
Qj(t)Gj_i iS b is calculated as 

fli =exp (— !— ] T f 0, ^— V 2<i<5, (55) 



7»-i 



7»-i 



whereas Qi 



is arbitrarily chosen. In this case, 

/AT-l \ JV-1 

= n ex P ^ - JJ r 1 0. 



i=l 



1 

7i 



(56) 



All comparative curves show an excellent match between 
the approximate solution and the Monte-Carlo simulation for 
the considered scenarios. 

The LCR curves, presented in Figs. 2, 4 and 6, manifest a 
behavior of a typical fading channel, since the average LCR 
increases with the threshold until some maximum, and then 
decreases. The maximized LCR and the maximizing threshold 
depend on number of hops and per hop SNRs 7$. The AFD, 
presented in Figs. 3, 5 and 7, also manifests a typical fading 
channel behavior, since AFD continuously increases with the 
threshold. 



Regardless on the per hop SNRs, for some specified thresh- 
old, the iV*Rayleigh fading signal typically remains less time 
in fading with the increase of N (Figs. 3, 5, and 7), whereas 
its LCR increases (Figs. 2, 4, and 6). These differences are 
more pronounced with the increase of the per hop SNRs. 
Similarly, for some specified threshold, the LCR increases 
with the per hop SNR, whereas its AFD decreases. Namely, 
for the considered range of values for the N and per hop SNR, 
$ increases with respect to both of those parameters, yielding 
to those observations. 

V. Conclusions 

In this paper, motivated by the fact that the multihop AF 
Rayleigh fading channel is a cascaded one, we present a novel 
and general analytical framework for the exact evaluation of 
important second order statistical parameters, as the LCR and 
the AFD of this channel. Moreover, we present simple and 
efficient closed-form approximations for the aforementioned 
parameters, using the multivariate Laplace approximation the- 
orem. The accuracy of the presented mathematical analysis 
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and the tightness of the proposed closed-form approximations, 
were shown by numerical examples and extensive Monte- 
Carlo simulations. The material presented in this paper can 
be efficiently used in determining the packet length selection, 
power and bandwidth allocation over multiple hops, and max- 
imum delay /latency requirement [12] [29], or determining the 
delay spread of frequency-selective multihop channels [30]. 
Finally, it could be useful in the study of the second order 
statistics of the cooperative diversity systems. 

Appendix I 
Laplace Approximation Theorem 

Let D be a possibly unbounded domain in the multidimen- 
sional space R n , u and h be real-valued multivariate functions 
of x = [xi, ...,x n ], and A is a real parameter. Consider the 
integral 



J(A) 



xeD 



u(x) exp(— Xh(x))dx. 



(1.1) 



If i) the integral J(A) converges absolutely for all A < Ao, 
ii) function h has an absolute minimum x = [x\, ■ ■ ■ ,x n ] at 
an interior point of D (this is turn implies that x is a critical 
point of h, i.e., V/i(x) = 0), and Hi) the Hessian matrix 

2 h 



dxidxj 

3 2 /i(x) 

dx{ 
d 2 h(5t) 



d 2 /t(x) 
d x\dx2 
d 2 h(Z) 
dx'i 



d 2 h{x) d 2 h(x) 



dx n dx\ dx n dx2 



dxidx n 
9 2 fe(x) 
dx2dx n 

d 2 h{yi) 



(1.2) 



is positive definite, then, for large A, 



J(A) 



n/2 



u(x) 



exp(-A/i(x)), (1.3) 



where det(-) represents the matrix determinant. The Laplace 
approximation theorem was originally proven by Hsu [23] for 
A — ► oo. It was observed in [24] that in many cases of interest 
the Laplace approximation performs very well even in sub- 
asymptotic cases where A remains small. 
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